"Polar" and "antiferromagnetic" order in / = 1 many-boson systems 
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In a system of interacting f — 1 bosons (in the subspace where the total spin in the z-direction 
is vanishing), we prove inequalities for the ground state expectation value of the density of spin- 
bosons. The inequalities imply that the ground state possesses "polar" or "antiferromagnetic" 
order when the quadratic Zeeman term q is large enough. In the low density limit, the inequalities 
establish the existence of a sharp transition at g = when q is varied. 
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Recent progress in cold atom experiments has led to a 
renewed interest in the low energy properties of interact- 
ing many-boson systems. Especially the systems of spin- 
ful bosons [ll-0| , in which hyperfine spin degrees of free- 
dom couple to many-body quantum physics of bosons, 
are expected to show a variety of nontrivial quantum 
magnetic phenomena. See Q for a recent review. It 
is a fascinating challenge in many-body physics to give 
a firm (hopefully rigorous) theoretical understanding of 
these phenomena [6|. The present work is a first step in 
such a direction. 

It has been argued since the early days 0, [1] that 
/ = 1 Bose-Einstein condensates with antiferromagnetic 
interaction have peculiar magnetic phases called "polar" 
and "antiferromagnetic" , in which global spin rotational 
symmetry is spontaneously broken Q. This symmetry 
breaking is of a "nematic" type (in the sense that the 
symmetry breaking in, e.g., the positive z-direction is 
the same as that in the negative z-direction). Interest- 
ingly the existence of nematic long-range order has been 
established rigorously in a class of 5* = 1 quantum spin 
systems ^ which can be regarded as limits of the / = 1 
Bose-Hubbard model [13, HI- The nematic phases in 
Bose-Einstein condensates and in spin systems should be 
related with each other. 

So far magnetic properties of interacting boson sys- 
tem are studied mainly by mean field theories, and few 
rigorous results have been obtained. Here we present rig- 
orous inequalities for the ground states (in the subspace 
where the total spin in the z-direction is vanishing) of 
fully interacting spinful boson systems which imply the 
existence of "polar" or "antiferromagnetic" order when 
the quadratic Zeeman term is large enough. In partic- 
ular we can show the existence of a sharp transition in 
the low density limit. We also discuss weak coupling the- 
ory, and place the present problem into the context of 
the standard picture of long-range order and symmetry 
breaking. 

We hope that the present study can become a start- 
ing point of further mathematical physical researches on 
quantum magnetic phenomena in interacting boson sys- 
tems. 



Definitions: We study a system of N bosons with 
(hyperfine) spin / = 1 in (or in a finite region of 
R'^), where N is even. Extension to bosons with higher 
spins is trivial. As usual we denote by ipaii") and ipl-{r) 
the annihilation and the creation operators, respectively, 
of a boson at r e R^ with (the z-component of) spin 
a e {-I-,— ,0}. These operators satisfy the commuta- 
tion relations [ip<j{r),'ip^^,{r')] — SaM'S{r — r') for any 
a, a' E {+,—,0} and r,r' G M.^. More rigorously, for 
any single-particle state (wave function) ip = (i^(r))^gR3, 
the operator '■= J (fir ip{r) ipl.{r) creates a bo- 

son in the state ip with spin a. We also let Pair) :— 
'0j-(r)V'CT(^) be the density of the spin a bosons, and 
Pir) = J2ae{+,-,o} Mr) the density. 

We consider the Hilbert space Hq, which is spanned by 
states of the form {HjLi i^tj [Vj]} ^vac, where $vac is the 
state with no bosons in the trap, ipj are arbitrary single- 
particle states (in a suitable space), and aj e {+,— ,0} 
are arbitrary except for the condition X^jLi '^j — 0- Note 
that we have restricted ourselves to the space where the 
z-component of the total spin is vanishing. This choice 
is realistic since our Hamiltonian ([T|) conserves X^jLi '^j- 
Experimentalists can prepare an initial state in Hq , which 
should evolve into a ground state in Ho- 

Let us write the full Hamiltonian as 

H = Ho + V-qNo. (1) 
The Hamiltonian of the non-interacting system is 

•' ™c7e{+-,o} 

(2) 

where U{r) is a suitable single-particle potential which 
describes the trap. 

We can consider an essentially arbitrary interaction V 
which is nonnegative and leaves the space invariant. 
To be concrete, however, let us take the most standard 
low energy effective interaction of dilute f — 1 bosons 
V = Vo + V2, where 

Vo-.^l f d?r d?r' vo{r ~ r')Al^{r, r')io,o(r, r') (3) 
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and 



V2 



'rd^r'v2{r^r') ^ Al^{r,r')A2.,n{r,r') 



m=-2 



(4) 

represent the interactions between two bosons with to- 
tal spin and 2, respectively. See, e.g., [5]. Here 
vo{s) and V2{s) are suitable nonnegative potentials, and 

AF,mir,r') Eo-,<t' (-^' '"k' f^') ^^('^)'0<t' (^') (where 

{F, m\a, a') denotes the Clebsh-Goldan coefficients) is the 
annihilation operator for the state in which the total spin 
and its z-component of two bosons (at r and r') are F 
and m, respectively 12 1. 

In the literature one often considers the local interac- 
tions with vf{s) = gpSis) where F = 0,2. Here, to 
avoid mathematical complications, we do not treat the 
delta function interactions. But we shall discuss the limit 
of short-range interaction after presenting the main in- 
equalities in Theorem 2. 

Finally —qNo is the quadratic Zeeman term generated 
by a uniform external magnetic field, where g S M and 
No — J d^r po{r) is the total number of spin-0 bosons. 
The standard Zeeman term is omitted since it has no 
effects within the space Hq as long as the magnetic field 
is uniform. 

Long-range order and spontaneous symmetry breaking 
in the weak coupling limit: To put the problem into the 
context of long-range order and spontaneous symmetry 
breaking in many-body systems, we shall discuss (almost 
trivial) weak coupling limit in some detail. 

Consider the system whose Hamiltonian consists only 
of the non- interacting part Hq. Since Hq does not flip 
the spin of each particle, the model has highly degenerate 
ground states within Hq (isj . 

Suppose that one adds to Hq a small spin-dependent 
interaction V which favors two particles to have vanishing 
total spin. This should be the case if uo('S) < V2{s). Then 
it is expected (and can be proved in related models IJ]) 
that the degeneracy is immediately lifted and the ground 
state in Hq is unique and has vanishing total spin. 

Let us consider this unique ground state in the limit 
where V vanishes, and denote it as ^q'^'^. Since we are 
now dealing with a non-interacting system, the A^-body 
ground state is written in terms of the creation opera- 
tors aj. :— J d^r ipQ{r)ip'^^{r), where ipo is the unique nor- 
malized ground state of the single-particle Schrodinger 
equation —{h?/{2m)}Aip{r) + U{r)ip{r) = e(p(r). It is 
also known that the information that the ground state 
has vanishing total spin is enough to uniquely determine 
<i>k°° as 



where C is a normalization constant IB" 12 1- This in 
particular implies that the probability Prob[iVo = n] that 
the number Nq of spin-0 bosons is equal to n behaves as 



Prob[7Vo = ' 
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2»{(n/2)!}2 



^- (6) 



Let us define the "order parameter" by O :— 
Nq — {N/3). The rotation invariance of $q°° implies 
($ftcc^O$frcc^ = 0, which also follows from ©. The 
peculiar power law behavior of ([6]) manifests itself in the 
fiuctuation of 6 as ~ (4/45)7V2, which 

becomes enormously large when N is large. 

Such a behavior of order parameter is indeed typical 
in a many-body state which exhibits long-range order 
but no spontaneous symmetry breaking [I8]. As a well- 
understood typical example, consider the S* = 1/2 an- 
tiferromagnetic Heisenberg model on the L x L x L cu- 
bic lattice with the Hamiltonian Hn = j') ' ^j' 
(where the sum is over all the neighboring sites), and 
define the Neel order parameter (or the staggered mag- 
netization) by On := ^^^(-l)-"^-"'^^''?']^^ (where we la- 
beled the sites as j = (ji, j2, js))- Although the present 
model exhibits Neel order at low temperatures and in 
the ground state, the ground state is unique for any 

finite L. The uniqueness implies ($qs' '^h'I'gs) = S^'^" 
ing no indication of the Neel order. The order manifests 
itself in the large fluctuation of Oh (or, equivalently, an- 
tiferromagnetic long-range order) as ($§g, {OH)^^Qg] 



> 



(const.) (L3)2 [19I - I21I 1. The ground state $gs exhibits 
long-range order but no symmetry breaking. To get a 
"healthy" ground state which exhibits both long-range 
order and symmetry breaking, one introduces the ficti- 
tious "staggered magnetic field" h. For ft, > let ^^o' 
be the ground state of Hn — hOu- Then it is proved [22[ 
that the spontaneous staggered magnetization satisfies 
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lim lim , Oh* 



H,/i\ 
GS / 



> 0. 



(7) 



One can say that an infinitesimally small symmetry 
breaking field h (in the infinite system) has triggered the 
symmetry breaking. 

The situation with our ground state ^q'^'^ is quite simi- 
lar [23I . Here the quadratic Zeeman term —qNo precisely 
plays the role of the symmetry breaking field. Indeed it 
is easily seen (for the Hamiltonian Hq — qIVq) that the 
ground state for any g > is the "polar" state 



<J,tlOC 

^pol • 
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free 
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(n:even) 



0(-2) 



-7l/2 



— I ^^^^ I 
2 ■ 2 • 



at)"(a;at_)(A^-n)/2$^ 



(5) 



which has (^^'■^f , A^o^poD = N, and the ground state for 
any (7 < is the "antiferromagnetic" state 0] 



iN/2) 



,(a^aL)^/^<l>v 



(9) 
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which has {¥lf,No<^ 



.^1?) = 0. The ground state © 
(in which O fluctuates widely) is in a sense unstable. 
Infinitesimally small q leads to stable states ([5]) or ([S]), 
which break the spin rotation symmetry explicitly. This 
is the reason that the mean-field theory predicts the polar 
or antiferromagnetic ground states even for g = 

We note in passing that the state which in the 
form of fragmented Bose-Einstein condensate i25] , is also 
rewritten in a "mean-field form" 



AF 



— (const.) Po{a\_ 



e2'«aL)^$v 



(10) 



where Pq is the projection onto Hq. Note that ([TU]) 
is independent of ^? e K. Here {al_ — e^"^^ a^J) j \f2 cre- 
ates the eigenstate of co-aB S'^^^ -I- sin 5'^^-' (i.e., the spin 
component in a direction orthogonal to the z-axis) with 
the eigenvalue 0. Since a suitable rotation in spin space 
brings aj into (a^,_ — e^'^aL)/'\/2, one might be tempted to 
conclude that the polar state ([8]) and the antiferromag- 
netic state (jni), (ITUl) are essentially the same state. But 
this is not the case since there is a projection in the 
two states are distinct {23]. 

We stress that the polar state ([5]) and the antiferro- 
magnetic state © appear as a consequence of a break- 
down of the rotational symmetry. The symmetry break- 
ing is triggered by the applied field —qNo, which breaks 
the symmetry in the z-direction. This direction is arbi- 
trary. If the symmetry is broken in the x-direction, for 
example, the polar and the antiferromagnetic states be- 
come (const. )(a^_ — at )^$vac and (const.){(a!'^ -I- aj -I- 

a^l)(a+ — a/S ao + ^-)}^^^^vac, respectively. We also note 
that the present symmetry breaking does not distinguish 
between, e.g., the positive and the negative z-direction. 
In this sense the present rotational symmetry breaking is 
nematic. 

Main results: Let us turn to interacting models and 
discuss our theorems. Although we believe that the 
ground state of the present model exhibits Bose-Einstein 
condensation, we do not make use of that property here. 
We shall focus on the ground state expectation value of 
po := / (fr po{r) = Nq/N, the ratio of the spin-0 bosons. 
We first make the following elementary observation. 

Theorem 1: Let <? = 0, and suppose that the ground 
state $GS of the full Hamiltonian ([T]) is unique and has 
vanishing total spin. Then we have ($gSj Po'J'gs) = 1/3. 

The unique ground state with vanishing total spin is 
expected for "antiferromagnetic" interactions character- 
ized by vq{s) < V2is). See {Tij. 

As before let (po{r) be the unique normalized ground 
state of the single-particle Schrodinger equation corre- 
sponding to ([1]). Let us define the effective volume V as 
the smallest constant which satisfies |93o(f')P < 1/^ for 
any r. We also let gp ■= J d^svpis) > for F = 0,2. 

Theorem 2: Let $gs be any ground state (in Hq) of 



the full Hamiltonian ([Ij. For any q > 0, one has 

go + 2(?2 N 



($GS,/5o$Gs) > 1 



6g V' 



and, for any q < 0, one has 

($GS,Po*Gs) < 



2go + g2 N 
12\q\ V 



(11) 



(12) 



The inequalities (fTTj) , (IT2]) imply that the ground state 
of ([T]) exhibits polar or antiferromagnetic order in the 
sense that ($gs, Po'&Gs) is close to 1 or 0, respectively, 
when \q\ ^ gN/V (assuming g c:^ g^ 52) |27|. Note 
that one cannot expect a sharp transition as in the free 
model discussed in the previous part. Since the ground 
state (within T-Lq) for a general nonvanishing V is ex- 
pected to be unique for any 5 £ M, the expectation value 
($GS, A'o'I'gs) should vary continuously with q when N 
is finite. 

To see a sharp transition (as q is varied), one should 
take a certain limit. The most challenging and interesting 
setup is to consider a model with suitable antiferromag- 
netic interaction and take the infinite volume limit (with 
a fixed density p) as in the case of the Heisenberg model 
(O. But this problem seems to be formidably difficult 
from the current technique of many-body systems. 

What one can conclude from the inequalities ([TT|) . (fT2|) 

is 



lim (^GSi/Oo^Gs) 

{N/V)iO 




(13) 



which establishes the existence of a sharp transition in 
the low density limit. Note that this result is valid no 
matter how strong the interaction is. This provides a par- 
tial rigorous justification of the prediction by the mean- 
field theory |5|]. 

Since the inequalities in Theorem 2 depend on the in- 
teraction only through go and .92, it is easy to consider 
the limit of short range interaction. One starts with po- 
tentials vq{s), V2{s) which are sufficiently "nice" so that 
the ground state $gs is well defined. After stating the 
inequalities ([TT1) . (fT2|) for ^gs, one may formally take the 
limit vpir) -> gFS{s) where = 0,2. The inequalities 
(ITT]) . remain valid as they are (provided that the 
limits of the left-hand sides are well defined) . 

Proof: To prove Theorem 1, we note that the condi- 
tions of the theorem implies that the ground state ^gs 
is invariant under rotation in spin space, and any rota- 
tionally invariant state $ satisfies (<I>,iVo$) = One 
way to see the latter fact is to rewrite $ in the language 
of distinguishable identical particles (by working out the 
symmetrization explicitly). Let Sj be the spin operator 
of the j-th particle. Since (Sj)^ — 2, the rotational sym- 
metry implies ($, (S'j^'')^$) = 2/3, which implies that 
the probability that the j-th particle has spin-0 is 1/3. 
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Theorems 2 is based on the following simple variational 
argument. For any G Ho, one has {^qs, H^qs) < 
(vj^ji/vj/). By using ([T]), we rewrite this as 



(?($GS,A^o4'GS> 



£;o-(*,i^o*)-(^',V'^'>, (14) 



where we noted ($gs,^^'gs) > and ($gs, -^o'I'gs) > 
i^o, where E'o is the ground state energy of Hq- 

For g > 0, we choose the variational state as the 
polar state $[j°f of dH). Then one readily sees that 

{'^',No^) = N and (\l/,iJo*) = Eq. The remaining 
(^f, can be evaluated explicitly as follows. Observe 
first that 

4o(r, r')io,o(r-, r') = Ul{r')i:l{r)Mr)Mr') + • • • , 



where the omitted terms involve spin + and 



(15) 
and has 



no contributions to (3>pof, K)*&pof)- By noting that 



pol 



Nl 



(16) 



we find 



(<fpoT,v^,<i>[,'oT: 



< 



i J d3rdVwo(r -/)(*',*') 

goN^ 



6V 



(17) 



where we used |(^o('"')P — 1/^ to get the final line. By 
evaluating ('I'pof , V2'&po°) m a similar manner, we get the 
desired inequality ([TT]) . 

For g < 0, we choose the variational state 4* as the 
antiferromagnetic state of ([9|), and proceed in the 
similar manner to get (fT2|). 

Discussions: We have proved the inequalities which 
shows that an interacting / = 1 boson system exhibits 
polar or "antiferromagnetic" order when sufficiently large 
quadratic Zeeman term is present. Although the conclu- 
sion may not be surprising, we stress that it has been 
established rigorously in a system with arbitrarily strong 
interaction. We hope that the present study becomes a 
starting point of further mathematical works on quantum 
magnetic phenomena in interacting boson systems. 

We have also placed the present problem of nematic 
order in boson systems into the larger context of long- 
range order and symmetry breaking in quantum many- 
body systems. That discussion as well as the mean-field 
analysis 0, 0, [1] suggests that the transitions to the po- 
lar or antiferromagnetic states are more drastic when the 



interaction has an antiferromagnetic nature [28|. This 
picture is consistent with the fact that the existence of 
nematic long-range order is proved in = 1 quantum 
spin systems which can be regarded as limits of the / = 1 
Bose-Hubbard model with antiferromagnetic interaction 
[§-[llj. Unfortunately our inequalities do not make any 
essential distinctions between antiferromagnetic and fer- 
romagnetic interactions. It is an interesting theoretical 
challenge to go beyond the present naive variational ar- 
gument to see essential roles played by antiferromagnetic 
interactions. 
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also thank Takuya Hirano, Hosho Katsura, Yuki Kawaguclii, 
Masahiro Takahashi, Tatsuyoshi Tanabe, Akinori Tanaka, 
Yuta Toga, and Masahito Ueda for valuable discussions. 
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